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ABSTRACT 



The conditions of masing and lasing in the hydrogen recombination Hnes (HLR) in the disk and 
outflow of MWC349 are studied. Comparison of the complete set of the observed a-lines, from U.2a 
through H92a, with simple models of optically thin spontaneous emission shows that observable HRL 
masing in this source is limited to the interval of the principal quantum numbers n « 10 36. 

We use our analytical and numerical results for the conditions of optimum masing (Strelnitski et al. 
1995b; "Paper 2" ), and the morphological parameters of a photoevaporized circumstellar disk modelled by 
Hollenbach etal. (1994), to obtain analytical approximations and quantitative estimates for the expected 
unsaturated maser gain and the degree of saturation in HRL masers in the disk of MWC349. It is shown 
that the unsaturated maser gains of the IR and optical HRL lasers should be very high but in fact all 
of these transitions are strongly saturated. Due to saturation, their relative intensity with respect to 
the spontaneous emission should steeply drop toward smaller n's which explains why optical and near 
infrared HRL lasers are not seen in MWC349. This result is quite general and makes the prognosis for 
observable high-frequency lasers in other sources pessimistic as well. An exception might arise due to 
special geometry, if lasing were confined to very small solid angles. 

We show that weak masing from the outflow of MWC349 is only possible in low n a-lines (n < 20- 
30), and then only if the outflow begins much closer to the central star than the Hollenbach etal. model 
predicts (« 3 • 10^^ cm). We therefore conclude that maser emission from the disk, rather than from the 
outflow, is responsible for the observed weak amplification of the "pedestal" spectral components in the 
a-lines with n<40. The fluxes in the lower frequency a-lines are well explained by spontaneous emission 
from the outflow, with proper corrections for free-free absorption. 

We review the current state (and prospects) of observations of masing and lasing hydrogen /3-lines 
in MWC349. We argue in favor of "close in n," rather than "close in frequency" a//3 pairs for intensity 
analysis, when maser amplification is present. 

Subject headings: stars:early-type — stars: individual (MWC349) — stars:mass loss — masers — radio 
recombination lines 
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1. INTRODUCTION 



MWC349, a well-known variable emission line star {V 1478 Cyg) with high IR excess and the bright- 
est stellar radio continuum source in our Galaxy, recently became the first known hydrogen maser in 
space. Martin-Pintado etal. (1989a, b) detected strong, double-peaked emission in its mm hydrogen re- 
combination lines (HRL), and even stronger submm double-peaked lines were subsequently detected by 
Thum etal. (1994a, b). Martin Pintado etal. (1989a) presented arguments for the maser nature of the 
detected lines and demonstrated the principal possibility of high-gain masing in these lines in the dense 
ionized outflow from the star, but failed to explain the observed spectral pattern kinematically. Pono- 
marev etal. (1989) attempted to explain this pattern, ascribing one of the spectral peaks to the outflow 
and another one to the partly ionized, infalling circumstellar disk. However, later, with new observational 
evidence, several authors (Gordon 1992; Planesas etal. 1992; Thum etal. 1992) were able to successfully 
ascribe both peaks to a circumstellar Keplerian disk. 

The masing Keplerian disk model has further been elaborated upon by Thum et al. (1994a, b) and by 
Ponomarev etal. (1994). The slight observed variations of the 35a-30a line widths with intensity, and the 
deep dip between the two spectral peaks of 30a, were respectively attributed by Thum et al. (1992) and by 
Ponomarev etal. (1994), to the effects of unsaturated laasev amplification. Earlier, however, Ponomarev 
etal. (1991) estimated that the 30a line should be at the threshold of saturation. Thum etal. (1992; 
1994a, b) found several observational manifestations of increasing saturation toward higher frequency 
submm lines and connected them with a decreasing role of coUisional relaxation with increasing line 
frequency. Saturation seems to be a key issue for predicting HRL lasing in the IR and optical; therefore 
it is one of the central topics of this paper and of the previous one (Strclnitski etal. 1995b; "Paper 2"). 

Originally the mm and submm HRLs in MWC349 were decomposed into a "masing" (double-peaked) 
component associated with the disk, and a "thermal" (pedestal) component attributed to the outflow 
(Thum etal. 1992). However, several recent papers claimed weak masing in the pedestal component too 
(Martin-Pintado etal. 1994a, b; Thum etal. 1994b; 1995). In the present paper we attempt to develop a 
self-consistent model of HRL masing and lasing in the disk and outflow of MWC349, in order to explain 
the full set of existing observations and to predict the properties of yet unobserved, potentially lasing 
HRLs at higher frequencies. 

We summarize, in Section 2, the observed fluxes in most of the hydrogen a-lines which have so far 
been detected in this source from Ha, at 6563 A, through 92a, at 3.6 cm. Comparison of the observed 
fluxes with simple theoretical models for optically thin spontaneous emission provides further evidence 
for masing, and delineates the limits of "active" frequencies. 

Section 3 is devoted to the interpretation of the double-peaked maser lines. Using the approach 
to saturation analysis developed in Paper 2, we obtain a reliable observational estimate of the degree 
of saturation for the 30a maser, conflrming the rough previous suggestions that this line should be at 
the threshold of saturation. We use the morphological parameters of the photoevaporation disk model 
of HoUenbach etal. (1994), and our analytical and numerical results concerning maser gains in HRLs, 
to calculate the expected values of the unsaturated gain and the degree of saturation for a-lines in the 
interval 5< n < 35. We show that all the IR a-line lasers in MWC349 should be strongly saturated, in 
accordance with the saturation tendencies revealed by Thum etal. (1994a, b) from their observations of 
mm and submm masers. The expected maser-to-spontaneous intensity ratios for these saturated lines 
steeply decrease toward lower n, providing further support for our previous conclusion (Strclnitski etal. 
1995a) that detectable double-peaked a-line lasing in MWC349 should extend from the submm into 
the IR domain — and then vanish somewhere around n w 10. We use the qualitative ideas about the 
interaction of saturated masing lines developed in Paper 2 to explain the observed "clustering" of the 
mm and submm lines near a single value of the separation of two spectral peaks in radial velocity. 

An analysis of the conditions for masing in the MWC349 outflow (Section 4) brings us to the 
conclusion that although weak masing in a-lines is possible here for n^30, the observed increase of the 
pedestal relative intensity toward lower n is more probably due to amplification in the disk. The bulk 
of the higher n pedestal lines is fairly well explained by optically thin, spontaneous emission from the 
outflow, when proper corrections for the free-free absorption in the lines arc included. 

In Section 5 we review the observational and theoretical situation of the masing hydrogen /3-lines. 
We demonstrate that, in accordance with our model predictions, the only known weakly masing double- 
peaked /3-line (32/3) forms in the disk near the corresponding a-lines of close n. We argue that, when 



3 



masing is present, observations of the "cilose in n" a./ [3 pairs should be more instructive than observations 
of "close in frequency" pairs. A prognosis for detecting other /?-hne masers in MWC349 is given. 

In Appendix 1, we argue that previous attempts to prove the presence of masing in mm HRLs in 
MWC349 are not fully satisfactory, and we give a more rigorous proof using two observed parameters of 
30a — the flux, and the interferometric angular scale of the emitting medium — and using only very 
general model considerations about maser geometries in a Keplerian disk. 

In Appendix 2, we derive a simple expression for the (3 / a line spontaneous emissivity ratio. It can 
be useful for estimating whether or not effects of optical depth (positive or negative) are present in the 
observed (3/ a intensity ratio. 

2. SUMMARY OF THE HYDROGEN a-LINE FLUXES IN MWC349 

In Figure 1 we plot the integrated intensities of the most of hydrogen recombination a-lines observed 
so far in MWC349 (we omitted some lines for the sake of the figure's clarity). 

The fluxes in clearly observed double-peaked, narrow-band components arc shown with dots; those 
in single-peaked, broad-band components ( "pedestals" ) — with squares. The latter symbol is also used 
for poorly resolved profiles and for profiles with only barely resolved double-peaked components. Filled 
symbols mark single observations (when several arc known, we normally choose the brightest one); open 
symbols, averages for some period of observations (taken from Thum etal. 1992; 1994a). Three optical 
and near IR lines (Ha, Pa and Ba) were corrected for interstellar and circumstellar extinction using the 
extinction curve given by Thompson etal. (1977). 

The two solid lines in Fig. 1 show the calculated spontaneous emissivity ratios for the Baker and 
Menzel's "Case B", for two extreme density values poossible in the emission line region of MWC349. 
The difference between the two lines is not significant. Since these theoretical lines show only ratios 
of emissivities, their position along the ordinate is arbitrary. The broken line gives the spontaneous 
emissivity ratios, corrected for free-free absorption, which is important for low-frequency (n > 30) lines. 
The correction procedure is described in Section 4. The position of this theoretical line along the ordinate 
is also arbitrary. 

Although the observational data in Fig. 1 are not homogeneous and belong to (or are averaged over) 
different epochs, wc can clearly sec some important tendencies. 

The most striking fact is that theoretical curve for optically thin spontaneous emission, with only the 
correction for free-free absorption, describes satisfactorily the observed fluxes in almost all the pedestal 
components, from the highest frequency optical and near IR lines to the lowest frequency radio lines — 
more than 4 decades in frequency and almost 12 decades in fiux! 

In contrast, the mm and submm double-peaked lines (dots in Fig. 1) lie, as a group, considerably 
higher than the spontaneous emission curve. This may be regarded as strong, though indirect, evidence of 
their being amplified. The pedestal components of the two submm lines, 26a and 21a, recently detected 
by Thum etal. (1994b), also show evidence of amplification in Fig. 1. 

To date, 10 double-peaked mm and submm hydrogen recombination lines have been reported in 
MWC349 (Marti'n-Pintado etal. 1989a,b; Thum etal. 1994a,b; 1995): nine a-lines — 36a, 35a, 34a, 
31a, 30a, 29a, 27a, 26a, 21a, and one /3-line — 32/? . The 36a line is the lowest frequency a-line where 
the double-peaked component first appears with intensity comparable to the pedestal component. Both 
the peak fiux density S of the double-peaked component and its integrated flux F increase with line 
frequency, attaining Smax ~ 300 Jy and f w 1 • 10"^° Wcm'^ by 21a (Thum etal. 1994b). 

The tantalizing question naturally arises: how far into the still unexplored IR domain does the 
phenomenon of double-peaked masing extend? There is one evident limit: the well known optical and 
near IR lines seem to be satisfactorily described by optically thin, or slightly thick, spontaneous emission 
(Fig. 1 and also Thompson etal. 1977; Hamann and Simon 1986). Smith etal. (1995) and Strelnitski etal. 
(1995a) recently undertook to narrow the gap of unobserved lines between submm and near IR. Smith 
etal. detected 6a (12.37 /xm) and 7a (19.06 yum) on the IRTF telescope, and Strelnitski etal. marginally 
detected 10a (52.53 /im) and obtained an upper limit for 15a (169.4 /im) with the KAO observatory. 
Fig. 1 shows that whereas the 7a line lies very close to the optically thin spontaneous curve, as other 
optical and near IR lines do, the 10a is w 4 times brighter than the spontaneous curve predicts. The 15a 
flux upper limit lies about 10 times higher than the spontaneous curve at this frequency. 
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Though Strchiitski etal. (1995a) cioukl not spcctrahy resolve 10a to see whether or not it is double- 
peaked, they conjecture, on the basis of the line's flux excess, that observable lasing from MWC349 might 
extend into the far IR down to n^i 10. However, at lower n detectable lasing definitely seems to vanish. 

An important characteristic of the masing mm and submm double-peaked lines is their relative 
intensity with respect to the pedestal. Along with the increase of the peaks' absolute intensity, the ratio 
Ic/Ipk of the pedestal's intensity at the central dip between the peaks to the peaks' intensity decreases 
from the threshold 36a line, where Ic/Ipk ~ 1, down to 30a, where this ratio falls to « 0.07 (Ponomarev 
etal. 1994). It is seen from the recent observations by Thum etal. (1994b) that Ic/Ipk is still low at 26a 
(w 0.06), but at 21a it becomes much higher, > 0.2. 

The following questions thus have to be answered: 

(1) Why does the double-peaked maser first appear at 36a? 

(2) Why does its relative intensity with respect to the pedestal component first increase and then 
decrease with decreasing n? 

(3) Why do the observable masers, both the double-peaked and the pedestal ones, eventually vanish 
with decreasing n? 

We address the double-peaked and the pedestal components separately in the two following sections. 



3. THE DOUBLE-PEAKED COMPONENT: DISK MASERS 

3.1. Morphological Model 

We analyse the conditions for HRL masing in MWC349 in the framework of the HoUenbach etal. 
(1994) morphological model of a circumstellar disk photoevaporized by a massive hot central star. Ac- 
cording to that model, in the case where the "fast" stellar wind from the hot central star has normal 
parameters, photoevaporation of the disk by Lc photons from the star creates two distinct regions: 

(1) the inner region, R < Rg = lO^^M^/M© (cm), — a nearly static ionized atmosphere of the disk 
whose scale hight H increases with the radial distance R as 

H{R) = Rg^-^y (3.1.1) 

and the base number density of electrons decreases with R as 

Neo{R) = 5.7 • 1029$^/2i?-3/2 cm-3 (3.1.2) 

in the bulk of the atmosphere, except in its innermost part where the density increases toward the center 
somewhat faster (oc i?~^/^); $49 here is the star's Lyman continuum photon luminosity, in 10*^ photon/s; 

(2) the outer region, R> Rg — the "disk wind," freely outflowing ionized gas. 

We assume for MWC349-A: = 26M0 (Ponomarev etal. 1994), thus i?g « 3 . 10^^ cm. The quantity 
$49 is more difficult to determinate for stars as peculiar as MWC349-A; therefore we normalize Eq. (3.1.2) 
on the 30a line. We do so by adopting 2R (30a) = 1.2 . 10^^ cm (the interferometric distance between the 
two hot spots, corresponding to the two spectral peaks, Planesas et al. 1992), and A^eO (30a) = 4. 10^ cm~^ 
(the optimum electron density for the amplification of this line at = IC^K; Paper 2). This fixes $49 
in Eq. (3.1.2) at the realistic value of $49 « 1.1. Thus, we adopt for the quasi-static part of the disk 
[dropping the sub-index "0," spreading thereby the Neo{R) value over the thickness H{R)]: 

Ne{R) = 5.9 • 10^9 cm-3 . (3.1.3) 

Combining this with the approximation for the optimum amplification density derived in Paper 2 (good 
over the whole interval 5,000 K< Tg < 10,000 K): 

iV™«^ « 8 . 10^^ (3.1.4) 
we obtain the following approximation for the radial distance of the maximum unsaturated amplification: 

i?max « 1.8 • 10^ n^-^ cm. (3.1.5) 
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The values of Rmax for the two extreme hues considercid here, 35q and 5a, differ by about 3 orders of 
magnitude (Fig. 2). They are, however, both located within the static atmosphere. 



3.2. Unsaturated Maser Gain 

In the interval of Rmax considered here the maximum amplification coherent path length (with 
respect to the radial velocity — the chord Imax in the Keplerian disk; cf. Paper 1) — is equal to Rmax 
to within a factor of order unity. Thus the maximum unsaturated optical depth can be estimated as 

max ^ imax i ^ imax p /'^ 1 ^ 

'0 ~ ^net '-max ~ l^net J^rnax ■ [O.Z,.l) 

Substituting into Eq. (3.2.1) Rmax from Eq. (3.1.5) and fc™r from Eqs. (3.1.4) and (3.1.5) of Paper 2, 
we find 

n -4.2 

^max _ _2 . iq2 _ ^j,^ ^ ^qA j^) ^ (3 2 2) 

To"''^ w -1 • 10^ ^"^'^ (Te = 5 • 10^ K) . (3.2.3) 

These equations show, first of all, that there should be a threshold principal quantum number (frequency) 
for a-line masing in the disk — the highest value of n = rithresh for which r™"^ first becomes ^ 1. 
Eq. (3.2.2) gives, for example, t™°^(36q;) « 0.9, and this is in good agreement with the observed threshold 
of the double-peaked masing at n « 36. The cut-off of the masing at n « 36 is, in this model, an "optical" 
effect based on the interplay of the density dependence of maser gain with the density structure of the 
disk; it does not require a physical cut-off of the disk as it does in Thum's etal., 1994a, suggestion. We 
note, however, that in fact 36a is masing in our model at ~ 1 • 10^^ cm — not so far from the edge of 
the static atmosphere {Rg « 3 ■ 10^'"' cm). 

Eqs. (3.2.2) and (3.2.3) show also that r™°^ steeply increases with decreasing n: from r™"^ ~ 1 
at n « 35, through Tq™''^ ~ 10^ at n « 20, through r^""^ - 10^ at n w 10, and up to r^"^ ~ 10^ at 
n w 5. Thus, if this photoevaporation disk model for MWC349 and the hydrogen population calculations 
are basically correct, the unsaturated maser gains in the IR and optical hydrogen recombination lines in 
MWC349 should be very high. 

Why th(ui are the observed fluxes in optical and near IR lines not extraordinarily high and are, 
in fact, close to the optically thin spontaneous theoretical curve, Fig. 1? To answer this question, we 
examine more closely the conditions of saturation for the disk masers. 

3.3. Saturation: Observational Evidence 

Ponomarev etal. (1991) briefly analysed the conditions of saturation in the H30a line in MWC349 
and showed that, with the observed flux and the geometrical scale of the emitting region estimated by 
Martin Pintado etal. (1989b) from temporal variations, this line is probably marginally saturated. 

A more detailed analysis of the connection between the observed line parameters and the degree 
of maser saturation was given in Paper 2. It was shown that the solid angle of maser emission — the 
most undetermined free parameter in the previous saturation analyses by Thum etal. (1992; 1994a) — 
can actually be excluded from final equations. Then the geometry of the problem can be reduced to 
some linear scale I, whose reliable estimate, for such a source as MWC349, can be taken directly from 
interferometric observations. The degree of saturation, as measured by the ratio of the actual average 
intensity in the source J to the "saturation intensity" Jg , was shown to be 



J. -O-'Wc) [To) (57)1— J • <^^^^'' 

where D is the distance to the source, S is the observed flux density, and I is the amplification path 
length. 

We apply now this condition of saturation to the double-peaked component of 30q! in MWC349, 
for which Planesas etal. (1992) measured the angular separation of the "blue" and "red" hot spots: 
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6s = 0.065 arcscc. It is convenient, in this case, to pass from / in Eq. (3.3.1) to Oi = (l/D) the angular 
measure of I as it would be seen from Earth, had it been perpendicular to the line of sight: 



T /\4.7/r.\/ \2 

J ^ ^ ^ -7 1 n \ I o \ I arcsec \ 



To a factor close to unity the amplification chord length in a Keplerian disk with our parameters is Z « s/2, 
thus, 6i w 0.50s ~ 0.03 arcsec. Substituting into Eq. (3.3.2) this value and also n = 30, S^oa ~ 25 Jy 
(Thum etal. 1994a), we find: J/Js ~ 0.5. This confirms the suggestion made before by Ponomarev etal. 
(1991) and by Thum etal. (1992; 1994a) that this line should be at the threshold of saturation. 

We emphasize that the forgoing analysis of 30q; saturation is based on the directly measured quan- 
tities — the flux density and the angular scale and is only weakly model dependent. Unfortunately, 
direct measurement of the angular scale has only been done for the 30a double-peaked component so far, 
and we can not apply the above analysis to other lines to see how the degree of saturation changes with 
n. However. Thum et al. (1992: 1994a, b) draw attention to some observational indications of growing 
saturation toward higher frequency lines: some broadening of the lines and levelling off of the maser 
photon luminosity. In the framework of their saturation analysis (which was reduced to the evaluation of 
the "saturation temperature" T, an analogue of J a) they find that maser saturation "becomes easier" 
toward lower n because Tg for them is smaller. These authors correctly notice, however, that such an 
analysis is insufficient: only comparison of a saturation parameter (Tg or Jg) with the corresponding 
expected parameters of radiation field in the source enables any statement about the degree of saturation 
in a masing line. 

In the theoretical analysis of saturation given in the next two sections we shall show that the degree of 
saturation indeed grows toward higher frequency HRL transitions in MWC349, and it leads to important 
consequences concerning observability of HRL lasers. 



3.4. Saturation: Theoretical Prognosis 

Assuming the solid angle of maser radiation Q, to be constant along the amplification path, and 
ignoring spontaneous emission with respect to the maser emission in a strong maser, the unsaturated 
growth of the average intensity with |ro| is given by 

Ju « ^^el-ol , (3.4.1) 

where /q is the input intensity — the intensity of the continuum or spontaneous emission, whatever is 
brighter, at the frequency of the line. In the Raleigh- Jeans domain (all the lines in question fall into 
this domain) the input intensity due to continuum is 2kT^ / < 2kTf. / A^, where is the continuum 
brightness temperature and Aq is the transition wavelength. It can also be shown (e.g.: Strelnitski 
1974) that the input intensity due to spontaneous emission is 2k\Tx\/ A^ where is the transition 
excitation temperature. Since |/3i2| ~ 1 at the density of maximum maser gain (Paper 2), we have 
\Tx\ ~ /|/3i2| ~ T'e- Considering that the exact value of /q is not very important when |r| is high, we 
can thus assume 

/o « ^ , (3.4.2) 
It was shown in Paper 2 that the saturation intensity for an HRL maser can be represented as 

where At = {A^-^ -|- A^^)~^ and Ct = {C^i -\- C:^)~^ are the harmonic mean values of the total Einstein 
A coefficient and of the total collision rate for the two maser levels, A21 is the Einstein A coefficient for 
the signal transition, and C21 is the rate of coUisional relaxation between the maser levels. It was also 
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argued that, at the optimum density for maser amplification in a HRL, At in Eq. (3.4.3) can always be 
ignored as compared with (Cj + C'21 ) which, in its turn, can be approximated as 

S = {Ct + C21) !v 3 ■ 10~^ Ne . (3.4.4) 

Substituting for A^e the approximation (3.1.4) for iV™°^; using: An+i.n ~ 6.3 • 10^ {^~^) z/q = 
c/Ao ~ 6.6 • 10^^ (s~^) (the common approximations valid for na-lines with n ^ 1); and neglecting 
At in Eq. (3.4.3), as compared with (Ct + C21), we obtain from Eqs. (3.4.1) - (3.4.4), at a line's optimum 
density: 

Js^^^^2.10-'n-'\ (3.4.5) 
c A21 

and 

^ «0.3n-i-3nel"°l . (3.4.6) 

•Js 

Note that Eq. (3.4.5) is valid over the whole temperature range considered here (5,000 K< Tg < 10,000 K), 

whereas the factor e'"^"' makes Eq. (3.4.6) rather strongly temperature dependent. 

The solid angle of maser radiation can be estimated in our disk model of MWC349 as 

O « "-^^"""^ « 0.5-^ « 0.9 • « 4 ■ 10-4 . (3.4.7) 



max 



It is assumed in Eq. (3.4.7) that the "vertical" size of the masing hot spot is « H, and its "horizontal" 
size is « 0.5Rmax- The latter seems to be a reasonable estimate considering the radial density scale in 
the model disk (limiting the ring of maximum amplification), and the gradient of velocities in a Keplerian 
disk, as is directly shown by the numerical modeling of maser amplification in such a disk (Paper 1). The 
last two equalities in Eq. (3.4.7) are obtained by substituting H from Eq. (3.1.1), with i?g w 3 • 10^^ cm, 
by assuming R = R,nax, and by substituting Eq. (3.1.4) for Rmax- 

Combining Eqs. (3.4.6) and (3.4.7), we get, for the case of MWC349: 



^ « 1-10-4 nO-^el^"! . (3.4.8) 

The exponential factor is decisive in Eq. (3.4.8). If the degree of saturation Ju/Js is indeed ~ 1 at n « 30, 

where also |to| ~ 6 (both these results seem to be reasonable based on comparisons of observations with 
our models), then Ju/Js should quickly become J> 1 at smaller n's, because of the rapid increase of |r| 
with decreasing n [cf. Eqs. (3.2.2-3)]. 

We therefore conclude that all the hydrogen IR and optical a-line lasers in MWC349 should be deeply 
saturated. 



3.5. Observability of IR and Optical Lasers 
It is shown (e.g.: Strelnitski 1974) that the output intensity of a saturated maser can be represented 

as 

Is-^Js\tq\, (3.5.1) 

where tq is the unsaturated optical depth. In the case of a homogeneous amplification path of length I, 
To for a HRL is given by 

T^ = {K + kf)l, (3.5.2) 

with k'^ being the continuum absorption coefficient at the line frequency and k!^ being the unsaturated 
line absorption coefficient. However, for those HRL's that have been observed as strong masers {n "z 35) 
fcg ^ ki (see Fig. 10 in Paper 2), hence we can ignore k'^ in Eq. (3.5.2) as compared with fcf and write: 

ro«fc?^ = ^^^iV{Mi2/, (3.5.3) 
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where Bi2 is the Einstein coefRcient. Av is the Une width, N{ is the LTE population of level 1, and 
bi, Pi2 are the usual departure coefficients. 

Since the masing "hot spots" in MWC349 are not spatially resolved as yet, it is more convinient to 
discuss their observed radiation in terms of the flux density rather than intensity, even in terms of the flux 
density integrated over frequency — for easier comparison with spontaneous emission. The integrated 
flux density due to the portion of the source having a surface area da and an outgoing intensity I at the 
line center, is 

dF = I da / D'^ , (3.5.4) 

where Av is the line width, and D is the distance to the source. 

Substituting Eq. (3.5.1) for /, Eq. (3.4.5) for Jg and Eq. (3.5.3) for tq into Eq. (3.5.4), considering 

that I da = dVs is the elementary volume from which the saturated maser emission comes to the observer, 
and integrating over the whole volume providing this emission, we get: 

Fs = JdFs = ^^ J N[hl3i2dVs . (3.5.5) 

In obtaining this equation we have also used the common relation between the Einstein's A and B 

coefficients. 

It is straightforward to obtain an analogous equation for the flux due to optically thin spontaneous 
emission: 

hp A21 



dF,p = — ^ / N^b2dVsp . (3.5.6) 



Deviding Eq. (3.5.5) over Eq. (3.5.6), we get: 

F, S J N[bil3i2dVs 

Fsp n A21 J N^b2dV,p 



(3.5.7) 



Of the three ratios composing the right-hand side of Eq. (3.5.7) two are obvious — the ratio of the 
integrated numbers of emitting particles and the ratio of the solid angles of emission (47r for spontaneous 
emission and for maser emission). The third ratio, S/A21 is the ratio of the rates of emission, per 
particle per second. The rate of the saturated maser emission is S because it is determined by the 
collisional "recycling" of the primary radiative-radiative pumping (cf Paper 2). If the observed maser 
amplification occurs at the optimum gas density, N'^T'"'^ , 5 is obtained by combining Eqs. (3.4.4) and 
(3.1.4): 

(5»2-10^n-°-^. (3.5.8) 
Using the approximation for A21 mentioned in Sect. 3.4, we find 

— Pi A ■ 10'^ n^-^ . (3.5.9) 
A21 

Eq. (3.5.9) shows that the factor 5/A21 strongly favors high n transitions. The physical reason 
for this can be roughly explained as follows. As is shown in Paper 2, N^°'^ is close to the density of 
thermalization, the density under which the net rate of collisional transitions across the levels in question, 
« 5 {k.Te/hvo), surpasses the net rate of spontaneous radiative transitions which is roughly ^ A21. Thus, 
close to the density of thermalization, 5/A21 ~ [kTe/hvo) oc rv^. This explains, very roughly, of course, 
why the 5/A2\ factor in Eq. (3.5.7) favors high n levels. 

The forgoing result is important for our prognosis of the observability of HRL lasers. If there is a 
gradient of density in the source, and hence the spontaneous emission in a line comes not only from the 
region favorable for lasing in this line, the ratio of the integrals in Eq. (3.5.7) can be <C 1. It is seen 
then from Eqs. (3.5.7) and (3.5.9) that at small n saturated laser emission can not significantly exceed 
spontaneous emission in principle, unless the solid angle of lasing is Q <C 1. This greatly narrows the 
conditions of laser detectibility: the source has either to be homogeneous in density (this density being 
just optimum for lasing in a chosen line), and/or the solid angle of the laser beam must be very small, 
requiring a special source geometry. Of course, to be detected in a small solid angle such a laser must also 
have a fortuitous orientation with respect to the observer. If these conditions are not fulfilled, the laser 
emission will either be lost in a bright spontaneous background or will miss the observer geometrically. 
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In the particular case of the model masing disk of MWC349, the solid angle of masing, given by 
Eq. (3.4.7), does decrease with decreasing n, but it decreases slower than S/A21 does [Eq. (3.5.9)], thus, 
after masers have become fully saturated, i.e. for n ^ 20, the contrast Fs/Fsp should drop with decreasing 
n, and at some value of n the laser component should vanish in the brighter background of spontaneous 
emission. We can crudely estimate this value of n as follows. Taking as Fg/Fsp for 21a the ratio of 
the observed maser flux, shown in Fig. 1, to the level of the spontaneous flux given in this figure by the 
solid theoretical curve, Fs/Fsp ~ 10, we find, from Fs/Fsp oc n^-^ [Eqs. (3.5.7), (3.4.7), and (3.5.9)], that 
Fs/Fsp drops to « 1 at n « 8. This agrees with the lack of a significant lasing component at n < 10 in 
Fig. 1. 

3.6. Saturation and the (Peak Separation) - (Line Frequency) Relation 

Proper accounting for saturation may help solve one more puzzle associated with the multi-line HRL 
maser in MWC349. It was observed that the peak separation dv increased with the decreasing n in 
double-peaked masing na-lines — from « 40kms~^ at 36a; and 35a to « 50kms~^ at 30a (Thum etal. 
1992; 1995). This seemed to be in accord with the idea that higher frequency lines should mase at higher 
density, thus closer to the center of the Keplerian disk where the rotation velocity is faster. However, 
attempts to explain the 5v{n) correlation quantitatively, assuming a regular density structure in the 
disk, met with difficulties. Ponomarev etal. (1994) argued that the observed ratio Sv{30a) /Sv{35a) 
would require an improbably steep radial gradient of density (a > 4 in A''e cx r~"), if each line mased at 
its density of maximum unsaturated gain. The recently discovered 26a and 21a (Thum etal. 1994a, b) 
are even closer in their peak separation to each other and to 30a (all have Sv « 50kms^"'^), so that very 
unacceptable density gradients (a ~ 10 — 70 ) are required to realize them at their densities of maximum 
unsaturated gain in a Keplerian disk. 

To explain the levelling off of Sv{n), Thum etal. (1994b) proposed that the disk rotation curve 
itself levels off inside some radius. No physical mechanism for such levelling off in a circumstellar disk 
is, however, proposed. Also improbable is the hypothesis that the ionized disk "finishes" abruptly at 
R ^ 40 A. U. from the star, where iVg ^ 10^ — lO^cm"'^, so that all the submm lines are forced to 
form at this high-density edge. At least, it is certainly not the case in the Hollenbach etal. (1994) 
photoevaporation disk model we adopt here. 

We hypothesize instead that the clustering of submm lines in some ring of the disk to produce 
close values of Sv, is due to their interaction via saturation. The forgoing analysis was based on the 
assumption that masers and lasers in different lines are formed independently from each other, at their 
densities of maximum unsaturated amplification. Yet, only for the lowest n lines are the regions of 
maximum unsaturated amplification well separated in the model disk of MWC349 — the higher the n, 
the more the amplification regions for adjacent lines should overlap (Fig. 2). We argued above that the 
disk masers and lasers with n $ 30 are affected by satiiration. A possible mechanism for a maser to change 
its location under the influence of another saturated maser was described in Paper 2. Saturation in a 
masing line increases the decay rate of its upper level, thereby facilitating energy sink for the adjacent line 
above. As a result, the latter is able to mase at a higher density than its density of maximum unsaturated 
masing, with the consequence that adjacent lines will cluster spatially. One can suppose that this effect in 
MWC349 is initiated by some lasing saturated IR transition below 21a, at a density lO^cm"^, where 
the overlapping of negative absorption coefficients for this group of transitions becomes large enough. A 
verification of this mechanism by direct numerical simulation is under way. 

4. THE PEDESTAL COMPONENT 
4.1. Masing in the Outflow? 

Some evidence for masing in the pedestal component of HRL, as opposed to the double-peaked 
component, has appeared lately. 

Martin-Pintado etal. (1994a) observed narrow, highly variable, and therefore presumably masing 
features in 39a and 35a at "non-double-peak" radial velocities which they ascribe to the outflow. In 
addition, Martin-Pintado etal. (1994b) argued that the observed flux in 66a (which, like all the other 
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lines with n ^ 36, has only a pedestal component) was c;onsidc;rably higher than the expected spontaneous 
emission, and hypothesized that it was due to maser amplification. 

Thum etal. (1994b) find that the pedestal in 21a is relatively stronger than in 26a and attributed 
this to effects of stimulated emission in the outflow. Their observation of these two lines is plotted in 
Fig. 1, and it is seen indeed that there is a significant excess of the flux, relative to the probable optically 
thin LTE value, for both the double-peaked and the pedestal components, especially in the case of 21a. 

Thum etal. (1995) find that the pedestal flux ratios of several (3 /a pairs of nearly the same frequency 
— from 26a/ 33(3 to 40a/48/3 — are considerably smaller than the expected optically thin LTE ratios, 
and explained the relative enhancement of the a-lines by the stimulated emission in the outflow. 

In view of this observational evidence, it is interesting to investigate the conditions of masing in the 
outflow of MWC349, as compared with its disk. The radio continuum data for this source are adequately 
explained by a model of spherically-symmetric expansion at constant velocity (Olnon 1975). Adopting 
the value of the outflow constant obtained by Cohen etal. (1985), Ne{R) R'^ = 0.81 • 10^^ cm"-'^ (which 
corresponds to the spherically symmetric, 50 kms~^ constant velocity mass loss of 1.2 • 10~^ Mq/jv), we 
will have for the radius R{Ne) at which the outflow has a given electron number density A^e: 

R{Ne) « 2.8 • 10^^iV-°-5 . (4.1.1) 

The values of A:™^'^ presented in Table 1 of Paper 2, were calculated for the thermal line broadening 
alone, at the temperature considered. To account for the line weakening due to the expansion of the 

outflow, these values should be decreased by a factor of « w/Avd, where w is the observed line width, 
due to expansion, and Avd is the thermal Doppler line width. Thus, in this case maser gain is 

l^mau^l ^ |^max| ^(jyrnax) ^Avd/w) . (4.1.2) 

If masing for each line occurs in its "optimum" layer {N^ = iV™°^), Eqs. (4.1.1) and (4.1.2), and 
Table 1 of Paper 2, with w/Avd ~ 4, yield the following estimates for the maser gain at the lowest 
conceivable temperature, Te = 5,000 K: |T|(21a) = 7; |T|(30a) = 1.8; |r|(40a) = 0.4; |T|(66a) = 0.02; 
|r|(92a) = 0.002. These values decrease with increasing Te: at Te = 10, 000 K, |r| is > 1 for only n < 21. 

These estimates show that maser gain > 1 in the outflow is only possible for relatively low n a-lines, 
n < 30 (±10), with the exact values depending on the temperature. 

At longer wavelengths (n > 30), Fig. 1 shows a gradually increasing depression of the pedestal fluxes, 
as compared with the "simple" optically thin LTE model. This model (solid lines in Fig. 1) assumes 
that radiation in all recombination lines comes from the same volume of space. Then the flux is simply 
proportional to spontaneous emissivity, which is oc n^^ for n ^ 1. However, as can be seen in Fig. 10 of 
Paper 2, the free- free opacity surpasses the line opacity at n 30, and their ratio then steeply increases 
with n (k'l in this figure should be reduced by the factor w/Avd « 4 — 6 to account for the velocity 
dispersion in the outflow). Therefore, for the long wavelength lines only emission from the volume located 
above the surface Tc ~ 1 will be observed. In the approximation of a spherically symmetric, constant 
velocity outflow with NqRI = 0.81 • 10^'^ cm~^ and Te = 10* K, the radius at which Tc = 1 is 

R{Te = 1) ^ i?o « 2.4 • 1016 (_^) -2/3 (cm) , (4.1.3) 

and the volume emission measure above the level Tc = 1 is 

dV = NoRl • 47r ^ « 8.0 ■ 10^*i?o ' « 3.4 • 10^8(^)2/3 ^ _ ^^ ^ ^^ 

Jro R ^j^z 

Thus, the luminosity and the flux in the lines with n ^ 30 should be proportional to n~^n~^ oc n~^. 
This prediction is shown in Fig. 1 with a broken line. It describes quite well the observed fluxes from the 
outflow for the whole range n « 30 — 92. The greatest deviation from this theoretical curve (a factor of w 4 
excess) is the observation of 66a by Martin-Pintado etal. (1994b). These authors concluded that maser 
amplification in 66a (of about this factor) takes place in the outflow and they support this conclusion by 
model calculations, assuming a mass outflow rate of 9 • 10~^Mq yr"^ 8 times higher than the value 
of Cohen etal. (1985) adopted here. With the latter, as argued above, the negative optical depth of the 
outflow in 66a is only ~ — 10~^, and masing with a gain > 1 is impossible. The 66a flux measured by 
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AltenhofF et al. is significantly lower and c;loHcr to our theoretical curve. More observations in this line 
are needed to clarify the possibility of time variations and masing. 

4-2. Masing Disk versus Masing Outflow 

We showed in the preceding section that only HRL with 30 ± 10 can have mascr gain 1 in the 
model spherically-symmetric outflow with the mass loss rate of MWC349. However, lines with n<40 
require A^e~5 • lO^cm"^ for their masing (see Table 1 in Paper 2). According to Eq. (4.1.1), these high 
densities would occur in the outflow only at i? ^ 1 • 10"'^^ cm. Yet, in the photoevaporation disk model 
adopted here, there is no outflow so close to the star: the quasi-static disk's atmosphere extends up to 
= i?g » 3 • 10^^ cm for a star of 26 Mg (Sect. 3.1). 

We conjecture therefore that the pedestal masing, especially masing with a gain > 1, as in the 26a 
and 21a lines, occurs in the disk of MWC349, rather than in its outflow. 

The optical thickness of the disk increases with decreasing n. It was demonstrated in Sect. 3.2 that 
a dominant double-peaked component first appears at ?a 36a where gain in the chords can be w 1. But 
when this happens, the lower n lines in the outer parts of the disk already have their gains approaching 
unity and can therefore undergo a weak amplification. For example, at the physical edge of the static 
atmosphere {R = Rg), where N,, « 4 • lO^cm""^, the 40a line would have r 1, if its amplification 
path length is ~ Rg-, its weak amplification could thus produce some structureless supplement to the 
spontaneous emission. Amplified emission from the different radial velocities within the disk (see Fig. 1 
in Paper 1) which form the "pedestal" component should progressively increase with increasing |t| and 
decreasing n. This is in accord with the decrease of the (/?/a)**"'"*/(/?/a)°''* ratio toward lower n, seen 
in the Thum etal. (1995) paper and suggestive of an increase in the a-line pedestal amplification with 
decreasing n. 

At still lower n, not only the chords (which give the double-peaked component) but also other parts 
of the disk acquire maser gain sufficient for high gain masing, and they will form an ever growing masing 
"pedestal" to the double-peaked component. This disk-origina,ted pedestal can be broad, because of 
the radial velocity dispersion in the disk. The contrast between the double-peaked component and the 
pedestal should decrease with increasing saturation. Our numerical modeling of the masing in an edge-on 
Kcplcrian ring (Paper 1) shows that in a fully saturated regime the center-to-peak intensity ratio may 
increase up to Ic/Ipk ~ 0.5. The increase of Ic/Ipk from Ki 0.06 by 26a up to ~ 0.2 by 21a, as observed 
by Thum etal. (1994b), is in good agreement with this model, given the evidence of a growing degree of 
saturation with decreasing n. We anticipate a further decrease of the pcak-to-ccnter intensity ratio with 
decreasing n , due to the increasing saturation in the disk, before spontaneous emission washes out both 
the double-peaked and the "pedestal" masing components, at n » 10. 

5. /3-LINES 

Five recombination /?- lines of hydrogen have been detected in MWC349 so far: 39/? (Gordon 1994), 
48(3, 38/3, 33/3, and 32/3 (Thum etal. 1995). Comparison of the intensities of a//3 pairs proved to be an 
effective tool in analyzing masing in a-lincs. Besides, /3-lines themselves show some evidence of increasing 
masing with increasing frequency, and it appears likely a real high-gain maser in the higher frequency /3- 
lines will be detected in the near future. These two aspects of studying /3-lines in MWC349 are discussed 
below. 

5.1. /3/a Intensity Ratio 

Gordon (1994) detected the first /3-line in MWC349 (39/3), and he was the first to use a /3/a intensity 
ratio (39/3/31a, the lines of approximately equal frequencies), as evidence of masing — in the 31a double- 
peaked component. 

Thum etal. (1994c) used the /3/a intensity ratios of several other equal-frequency pairs they had 
detected, and demonstrated weak masing in the a-line pedestals (see Sect. 4). 
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The use of the (3/a ratio is based on the assumption that both hnes of the pair are produced in 
approximately the same region of space, and on the fact that \k^^^\ is normally several times higher than 
fcfgj. Because of the latter condition, the medium always becomes optically thick in the a-line first. 
Note that the relative effects of optical thickness are just opposite when the populations are inverted and 
when they are not. In the former case the (i/a intensity ratio changes in favor of the a-line, because of the 
higher maser amplification. In the latter case the ratio changes in favor of the /?-line, as compared with 
the optically thin case, because of the higher self-absorption in the a-line. Since in the /3/a method the 
observed intensity ratio is compared with the spontaneous emissivity ratio, it is useful to have a simple 
approximation for the latter. We derive such an approximation in Appendix 2. 

For the two pedestal pairs, 39/3/40a and 39/3/31q!, Gordon (1994) measured intensity ratios of w 0.6 
and 0.12, respectively. Using, as the corresponding spontaneous emissivity ratios, e(39/?)/e(40a) « 0.6 
and e(39/3)/e(31a) « 0.3, he came to the conclusion that these lines were observed in an approximately 
LTE optically thin ratio. Our Eq. (A2.5) gives for the e(39/3)/e(40a) « 1.3 — twice as much as the 
value adopted by Gordon. Comparison of our emissivity ratios with the flux ratios observed by Gordon 
leads us to the conclusion that both a-line pedestals are amplified by a small factor of « 2. We note 
a remarkably good agreement of our (/3/a)'*^'°"*/(/3/a)°^* ratio for the 39/3/40a pair with the Thum's 
etal. (1995) ratio for the 48/3/40a pair: w 0.5 in both cases. This is another confirmation that the 40a 
pedestal is amplified by a factor of w 2 (see Sect. 4.2). 



5.2. (3 -Line Masing in the Disk 



The highest frequency /3-line detected so far in MWC349, 32/3, is also the only one demonstrating a 
double-peaked component along with the broad pedestal component (Thum etal. 1995). By analogy with 
the double-peaked a-lines, the double-peaked component of 32/3 was interpreted by Thum etal. (1995) 
as a result of weak masing in the disk. Here we develop this interpretation further. 

First of all, using the rather precise information about the peaks' radial velocities (Table 2 in Thum 
etal.), we argue that the 32/3 peakes were formed in the same region where the a-line with the closest 
observed n (34a) was formed, rather than in the region where the 32/3 frequency companion (26a) was 
formed. The red peak of 32/3 (33.6 ± 2.3kms~^) is closer to that of 34a (32.3 ± l.Okms"^, as averaged 
for a long period of observations; Thum etal. 1994a) than to that of 26a (30.7 ± .1 kms~^), though the 
difference (2-3 kms""'^) is comparable with the observational errors or time variations (1-2 kms^^). Yet, 
the blue peak of 32/3 (—13.3 ± 1.5kms~^) is much closer to that of 34a (—13.0 ± 1.5kms~^) than to 
that of 26a (— 21.1 ± .lkms~^), and in this case the difference between the two groups (w 6kms~^) is 
significantly greater than the uncertainties (±2kms^^). 

If the two lines are indeedcoming from the same volume of space, masing seems to be necessary to 
explain their observed flux ratio. From Eq. (A2.5) we have: 

^»1.6. (5.2.1) 

e34a 

Yet, the fluxes in the peaks of 32/3 (2 ± 2 Jykms"^) and 34a (20 ± 20 Jykms^^) have a ratio of only 
--^ 0.1, though with a high uncertainty. The theoretical optically thin ratio (5.2.1) is thus much higher 
than the observed ratio; this discrepancy is naturally explained by maser ampliflcation which is stronger 
in 34a than in 32/3. 

That masing a- and /3-lines with close n can indeed share the same region of space, can be shown 
in the following general manner. Every three adjacent levels, 1, 2, and 3 (Fig. 3) produce two a-lines 
(2 — 1 and 3 — 2) and one /3-line (3 — 2). The definition of excitation temperature [Eq. (2.2.1) in Paper 2] 
and the obvious identity {Nig2/N2gi) (A^253/-^3S'2) {Nsgi/Nigs) = {ni/n^} {n2lnz) (na/ni) = 1 give the 
following relation between the frequencies and excitation temperatures of the three transitions (Strelnitski 
1983): 

7f — = Tf, — + Tf, — • (5-2.2) 

-1x31 J-x2\ -tx32 

Since /3-coefRcients of the adjacent a-lines, and, therefore, the excitation temperatures T^2i and T'x32) are 
always close to each other [see eqs. (2.2.3) and (2.2.4) and Fig. 4 in Paper 2], and since i/31 = 1^21 -h ^32, 
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it follows from Eq. (5.2.2) that 

Tr rri ^31 

x31 = J-X21 



V2l + 1^32 y". 



should also be close to Tx2i and TxS2- Thus, the (3{Ne) and the ki (N^) profiles should be similar for the 
a- and /3-lincs of close n. The knet (-^e) profiles will also be similar for those lines with ki fee (n S 40; 
see Fig. 8 in Paper 2). 

Fig. 6 in Paper 2 shows indeed that the knet (-^e) profiles are quite similar for the a- and 0-lmes of 
close n. The k^Jl/k^J^ ratio, at the density of their maximum gain {Ne ~ 1.5 • 10''cm~^) is w 0.2, at 
Te =5,000 - 10,000 K. Suppose, r^^f « -3. Then r^^f k. -0.6, and their intensity ratio, produced by 
unsaturated amplification, will be ~ e'^'^ w 10, as observed. According to Fig. 5, Paper 2, to produce 
these optical depths at 10, 000 >Tg> 5,000 K an amplification length of 50 100 A.U. would be sufficient, 
again in accordance with the interferometric distance of 80 A.U. between the 30a hot spots (Planesas 
etal. 1992). 

If the model of a masing Keplerian disk, with density increasing toward the center, is correct, then 
32/3 should be the "threshold" of masing for /3-lines, as 36a is for a-lines (Thum etal. 1995). The value 
of the maximum \knet \ is about an order of magnitude smaller for a /3-line than for an a-line of close n, 
which is why at the radial distance where the "threshold" a-line is formed the corresponding /3-lines are 
still optically thin and thus not amplified enough. The threshold /3-line, with |T„et| ~ 1, should appear 
at a smaller radial distance (higher density), where the corresponding a-line is already optically thick. 
This explains why the threshold of double-peaked masing occurs at smaller n for /3-lines (32/3) than for 
a-lines (36a). 

The intensity of double-peaked masing in /3-lines should increase further with decreasing n, as hap- 
pens among a-lincs. Using Fig. Qg of Paper 2 and Eq. (3.1.5) of this paper, we find, for example, that the 
maser gains of the 29/3 and 26/3 lines (both having good conditions for ground-based observations) should 
be, respectively, ~ 1.5 and 2.5 times higher than the gain in 32/3. Hence, maser amplification several 
times stronger in 29/3 and 26/3 lines is anticipated. Detection of a masing double-peaked component in 
these lines, which have very close observed a-companions (29a — 30a and 26a — 27a, respectively), will 
be especially useful for a comparative analysis similar to that presented for the 32/3/34a pair in Thum 
etal. (1995) and in this section. 

Suborbital and orbital IR observatories have a good chance of detecting, in the near future, even 
stronger lasers in higher frequency /3-lines — from w 20/3 through w 15/3 (w 200 — lOO/im), where the 
unsaturated maser gains can be considerably higher than unity. 

In conclusion, we emphasize once more that an important difference exists between the traditional 
analysis of optically thin a//3 pairs and the new situation where masing is responsible for the observed 
intensities: a "close frequency" pair is preferable in the first case, while a "close n" pair — in the second. 



6. CONCLUSION 

MWC349 is the only source in which strong HRL masers have been detected; moreover, it is a stellar 
source in which most HRL a-transitions (both masing and non-masing) have been observed more than 
twenty! This provides us a unique opportunity to study in detail the conditions needed for high gain 
HRL masing in the larger general context of HRL emission. Such a study was undertaken in this paper 
and Paper 2, and can be summarized as follows: 

(1) The factor iV| plays a decisive role in determining the gain |r| of an HRL maser. 

(2) At high principal quantum numbers n, the practical creation of a high gain HRL maser is hindered 
by the large dimensions needed to make |t| > 1, a requirement caused by the low densities involved in 
inverting high n transitions, and by free-free absorption which reduces the gain. 

(3) Low n transitions require high densities for inversion, and if present in a source (as is the case for 
MWC349), the unsaturated gain may become very large even with small dimensions. However, masers 
or lasers at low n saturate easily, and tend to disappear in the background of the spontaneous emission 
of the source unless the solid angle of the maser beam is very small. 

(4) In the case of MWC349 the optimum interval for strong a-line masing and lasing lies between 
these two extremes, 10 t 36. Relatively strong, double-peaked IR a-line lasers may, therefore, be 
discovered in this source in the interval n w 10 — 20. Their expected fluxes can be estimated from Fig. 1. 
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(5) A |r| > 1 Q-linc mascr in the outfiow of MWC349 doesn't seem possible for n Z 30 40. 

(6) Fig. 1 confirms the hypothesis of Thum etal. (1994b) that the relative increase of the pedestal 
component in 26a and 21a is due to masing. But we show the observed pedestal masing probably occurs 
in the disk, not in the outflow. We predict further increase of the pcdestal-to-pcak flux ratio toward lower 
n, until at n « 10 the whole structure is lost in the steeply growing spontaneous emission of the source. 

(7) High-gain, double-peaked /3-line masers and lasers with n < 32 should be present in the submm 
and IR domain, with 29/3 (614 yum) and 26/3 (447.4 /um) being the best candidates for the ground-based 
searches. 

(8) (3 /a line pairs of close n transitions are more informative in the case of masing than the close- 
frequency line pairs. 
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APPENDIX 1: Prof of the Need for Masing 



Though the hypothesis of maser amplification in the mm and submm HRL in MWC349 is now 

widely accepted, paradoxically enough no strict proof for this has yet been given. In this Appendix we 
briefly review the existing attempts to prove maser action occurs in these lines. We demonstrate their 
insufficiency, and propose a more rigorous proof. 

The high intensities of the double-peaked lines — higher than one would expect if the emission 
were purely spontaneous — was the primary reason to suspect maser amplification (Martin-Pintado 
etal. 1989a). These authors demonstrated, by model calculations, that the observed intensity of 30a, 
for example, could not be produced by thermalized emission from the expanding envelope of MWC349, 
the location where the lower frequency recombination radio lines were thought to arise. However, it is 
known now that the double-peaked lines do not originate in this outflow, so this modeling can no longer 
be taken as a proof of their masing. 

Martin-Pintado etal. (1989b) made an attempt to guess the size of the emitting "hot spots", with 
the (rather arbitrary) assumption that the upper limit of the size is deflned by the observed time scale 
of the intensity variations and the maximum observed velocity of the gas motions around the star. They 
obtained a size of :S 4 • 10^^ cm which, at the distance of 1.2 kpc (Cohen etal. 1985), corresponds to 
an angular size of < 0.022 arcsec. The observed flux densities then give a brightness temperature in the 
hot spot of Tfe > 1 • 10^ K implying some amplification, because the highest conceivable brightness 
temperature of the input radiation for the maser is less then 10^ — 10° K. 

Two more arguments were put forward in support of masing in these double-peaked lines: their 
time variability (Martin Pintado etal. 1989b; Thum etal. 1992), and the narrowness of their spectral 
profile (Thum etal. 1992; Ponomarev etal. 1992). Both phenomena are more easily explained with the 
hypothesis of unsaturated maser amplification. 

The last two arguments are, however, indirect, while the first two are based, respectively, on the 
incorrect association of these lines with the outfiowing gas, and on a mere guess of the radiating region's 
size. Thus none of them can be taken as a strict proof of masing. Yet, the direct measurement of 
the angular separation between the two hot spots responsible for two spectral peaks in 30q! (Planesas, 
Martin-Pintado, and Serabyn 1992) actually enable us to prove masing quite rigorously, as follows. 

According to Planesas etal. (1992) the angular distance between the centers of the two hot spots is 
9s = (0.065 ± 0.005) arcsec = (3.2 ± 0.2) 10~^rad. This separation can serve as a strict upper limit to 
an individual hot spot's angular diameter, 9d. If the separation between the two spectral peaks is due to 
the motions of the gas (which seems very probable), there should be some connection between the radial 
velocity scales and the geometrical scales. In particular, the narrowness of the peak's spectral width. Aw, 
as compared with the separation between the peaks, dv, {Av : « 1 : 4), is an indication that the ratio 
of the linear diameter of a hot spot to the distance between the hot spots, d : s = 6d '■ Og, should also be 
about 1 : 4. Our numerical modeling of maser amplification in an edge-on Keplerian disk (Ponomarev 
etal. 1994) indeed gives d : s ^ 0.3. Thus, we can safely admit that at least d^ <C (or 0^ <i: 0^) is true. 

Suppose that the observed flux F within one peak of the 30a double-peaked profile were due to 
isotropic, optically thin spontaneous emission. In this case 

where D is the distance to the source; 

L,p = N2A2ihuoV (A1.2) 

is the total spontaneous luminosity of the medium volume V with the upper state population N2; ^21 

is the Einstein coefficient for spontaneous emission, i>q is the frequency of the transition, and h is the 
Planck constant. Unless the geometry of a hot spot is very unusual, its "vertical" dimension should also 
be w d (or ^d, as it can be in the case of an edge-on disk). Then 

y < < s2^, or F < {OsOfl , {Al.3) 
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where / is the extension of the emitting region along the hne of sight. Combining equations (Al.l) - 
(A1.3) leads to the following strict limit for the integrated flux due to spontaneous emission: 

F -^^elNjlhvo . {AlA) 
47r 

The population density of a hydrogen atom level with principal quantum number n is given by the 
Saha-Boltzmann equation: 



Nn=(,r^] N^e^bn, (^1.5) 



^2 \ 3/2 

\ 2'KmkTe 

where N^. is the electron number density; /„ is the ionization energy of the state n; 6„ is the departure 
coefficient; k is the Boltzmann constant, and m is the electron mass. For n « 30, 6„ is close to one if 
Ne Z lO^cm""^, and for the temperatures considered here, 3000K < Te < 10000 K the factor e^ is very 
close to one too. Thus, we have for the n = 30 level: 

iV„ « 4 • IQ-^^N^ . (^1.6) 

Substituting into Eq. (A1.4): ^31^30 « 225 s"-'^, together with the Eq. (A1.6) and with the observed 
value Faoct « 1 • 10~^'*^;^§^ (Thum etal. 1992), we obtain finally the following strong requirement for the 
emission measure: 

N^l > 1 • 10^1 cm-'"^ . {A1.7) 

For a Keplerian disk l^s (Ponomarev etal. 1994). With the above intcrferomctric value of 6s and with 
the distance to the source D = 1.2 kpc, we have l^s = 6gD = 1.2 • 10^^ cm. Then another strong 
requirement, for the electron density, follows from Eq. (A1.7): 

ATe > 9- lO^cm"^ . {A1.8) 

We emphasize that (A1.7) and (A1.8) are strong inequalities, becose of the probable small size of d, as 
compared with s. But even these lower limits for N^l and A^e are high enough to show the model of 
optically thin spontaneous emission is internally contradictory. 

Indeed, if the 30a line were due to optically thin spontaneous emission alone, the observed line width, 
Av w 12kms,~^ corresponds to a kinetic temperature T,, « 3000 K. The Rayleigh-Jeans approximation 
is valid for mm lines, and the optical depth in a line of spectral width Ap is given by (Lang 1974): 



7re^ ^ Ni hu 
mc Av kTx 



where is the excitation temperature of the transition. At the densities set by Eq.(A1.8) the 30a 
transition should be deeply thermalized (see Fig. 66 in Paper 2). With Au » i/q ■ ^ » 9 • 10^ Hz (and, 
correspondingly, T^c = ~ 3000 K), and with the oscillator strength for hydrogen a-transitions being 
0.194n (Menzel 1969), we have 

Tij « 8 • IQ-^'^N^l . {AlAO) 

Substituting (A1.7) for N^l gives r » 1 which contradicts our assumption of optical thinness; thus this 

explanation of the observed value of F (30a) is contradicted. 

Let us suppose instead that the 30a line is thermalized and optically thick. In this case the flux 
density in the line center is 

Sma. « B{uo, T,) ■ ej < B{uo,T,) ■ el , (^1.11) 

where B{iyo,Ta;) is the Planck function. Substituting into Eq. (Al.ll): B(232 GHz, 3000 K) « 5 ■ 
10~^^ erg/(cm^sHzsr), and 0l « 1.0 • 10~^^ sr, we obtain, as an upper limit for the possible flux 
density: Smax 0.5 Jy, much less than the observed value, 5 w 25 Jy. 

We conclude that neither optically thin nor optically thick spontaneous emission alone can explain 
the observed flux density and integrated flux in 30a, and, therefore, maser amplification of the line in an 
optically thick inverted medium is necessary. 
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APPENDIX 2: Spontaneous Emissivity Ratios for /3/a Line Pairs 

The integrated spontaneous emissivity in the hne, whose lower and upper states' principal quantum 
numbers are, respectively, I and u, is 

eui = ^NuA^i , (^2.1) 
where h is the Planck constant, uq is the transition frequency, Nu is the upper state level population, 

^„i = 7.5-10-22i/oV;«(^) {A2.2) 

is the Einstein coefRcient, and /j„ is the oscillator strength. 

The transition frequency is given approximately by 2cR/l^ for a-transitions and AcR/l^ for /?- 
transitions, with R being the Rydberg constant (Lang 1974). Thus 

^o/3~2fM'. (^2.3) 



^13 

We emphasize that the values of I and u for the a transition may differ from those of the /3-transition in 
an a/p pair. 

When the upper levels of the a- and /3-transitions are different, but not too distant, we obtain from 
Eq. (Al.5), ignoring the small difference in the factor 6*=^= : 



With flu « 0.194Z for a-hnes, « 0.027/ for /?-lines (Menzel, 1969), we have finally from Eqs. (A2.1) 
through (A2.4): 

!^«l.lfM'^. (A2.5) 

If both la and l/s are ^ 1 and \la — Ip] i 10, the uncertainty of this approximation is ^ 10%, even if the 
bn ratio is set equal to unity. 

For the particular case of a pair of a- and /3-lines of close frequencies, the condition UQa ~ i^op, 
combined with Eq. (A2.3), implies l/s/la ~ 2^/^, and we get from Eq. (A2.5): 

^ « 0.3 « t^oa) ■ {A2.6) 
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FIGURE CAPTIONS 



Fig. 1. Observed fluxes in hydrogen recombination a-lines in MWC349. Dots: the double-peaked 

component. Squares: the pedestal component, or an unresolved line. Filled symbols — single measure- 
ments, either the only one-two known, the most reliable, or the brightest in a series. Open symbols — 
averages for some period of time (Thum etal. 1992; 1994a). The observed values of flux for 2a, 3a, and 
4a were corrected for extinction, in accordance with the Fig. 5 of Thompson etal. (1977). Error bars re- 
flecting internal convergency are not shown: they are comparable to the symbol sizes for all observations. 
"Error bars" by Ha show the amplitude of variations observed by Greenstein (1973) during one day (they 
were not intrinsic variations in MWC349-A — see footnote ^ in Brugel and Wallerstein, 1979). Two solid 
lines give, at two extreme values of N^, the flux ratios predicted by the optically thin LTE model ("Case 
B", Te = 10,000 K). The dashed line takes into account free- free absorption at low frequencies. Vertical 
position of all theoretical lines is arbitrary. 

Fig. 2. Location, within the model disk, of the annuli of maximum unsaturated amplification for na 

lines. The values of n are indicated along the top "border" of the disk. The two "borders" give, in fact, 
only an idea of the disk's flare, showing the Log of the scale hight (in 1.8 • 10^" cm) as a function of radius 
(in 1 • 10^^ cm). One quarter of the disk is only shown. The width of every annulus (at half maximum 
of fc™^^) is w ±0.2 in logarithm. Hence, in this logarithmic scale all the annuli have approximately the 
same width, and it is seen that overlapping of annuli should increase with increasing n. 

Fig. 3. To the explanation of the closeness of excitation temperatures (inversions) in a- and /3-lines 
of close n. 
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